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NOMINCLATURH 
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AR aspect ratio,/) 2 A 
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pressure gradient parameter. 
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boundary-layer displacement thickness, 
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molecular viscosity 
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H uid density 
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Subscripts: 

e edge of boundary layer 
f flap 

i inner region of boundary layer 

m main airfoil 

o outer region of boundary layer 

s slat 


free stream 


LEADING-EDGE SLAT OPTIMIZATION FOR MAXIMUM AIRFOIL LIFT 


Lawrence E. Olson, Phillip R. McGowan.* anti Clayton J. Guest* 
Ames Research Center 

SUMMARY 


A numerical procedure lor determining the position (horizontal location, vertical location, 
and deflection) of a leading-edge slat that maximizes the lift of multielement airtoils is presented. 
The structure of the flow field is calculated by iteratively coupling potential How and boundary- 
layer analysis. This aerodynamic calculation is combined with a constrained function minimization 
analysis to determine the. positjor.. of a leading-edge slat so that the suction peak on the nose of the 
main airfoil is minimized. The slat position is constrained by the numerical procedure to ensure an 
3 » tied boundary layer on the upper surface of the slat and to ensure negligible interaction 
between the slat wake and the boundary layer on the upper surface of the main airfoil. The highest 
angle of attack at which this optimized slat position can maintain attached flow on the main airfoil 
defines the optimum slat position for maximum-lift. The design method is demonstrated for an air- 
foil equipped with a leading-edge slat and a trailing-edge, single-slotted flap. The theoretical results 
are compared with experimental data, obtained in the Ames 40- by 80-Foot Wind Tunnel, to verify 
experimentally the predicted slat position for maximum lift. The experimentally optimized slat 
position is in good agreement with the theoretical prediction, indicating that the theoretical proce- 
dure is a feasible design method. 


I. INTRODUCTION 


The multielement wing is essential for generating the high lift required of modem fixed-wing 
aircraft during takeoff and landing. The necessity to carry high lift results in flow fields that are 
strongly dependent on viscous effects. In the past, the design of these configurations has been 
dependent on expensive and time-consuming wind-tunnel testing to optimize the placement of the 
various wing elements of a high-lift system. This task is formidable because each of the three param- 
eters (horizontal location, vertical location, and deflection) which specify the position of each flap 
and slat relative to the main airfoil must be optimized to yield the best aerodynamic performance. 
If this process could be augmented by theoretical analysis, the saving in time and wind-tunnel 
testing would significantly enhance the versatility and degree of optimization achievable in the 
development of mechanical high-lift systems. 

Stevens et al. (ref. 1), Bhateley and McWliirter (ref. 2), Callaghan and Beatty (ref. 3), and 
Olson and Dvorak (ref. 4) present theoretical methods for the analysis of multielement airfoils. 
Although these analyses (all of which rely on coupled boundary-layer and potential flow analyses) 
differ somewhat in the methods used to determine the structure of the viscous flow, references 1 
through 4 have demonstrated that, for flows in which the effects of boundary-layer separation are 
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negligible, it is possible to predict the lift of multicomponent configurations with reasonable accu- 
racy. These successes provided the impetus for utilization of these types of analysis techniques for 
the theoretical optimization of high-lift configurations. This optimization capability would be used 
to augment and, in some cases, replace experimental optimization. 


The problem of designing multicomponent airfoils has been addressed by several researchers, 
and the following arc several of the more significant examples. A design technique for shaping 
leading-edge slats to give a specified pressure distribution on the nose of an elliptical airfoil was 
developed by O'Pray (ref. 5). In his analysis. O’Pray requires that the slat position remain fixed and 
assumes that viscous effects are negligible. James (ref. 6) and Kennedy and Marsden (ref. 7) also 
present inviscid methods tor generating multielement airfoils with prescribed pressure distributions, 
but their methods do not address the optimization problem. Onnsbee and Chen (ref. 8) have con- 
sidered the shaping ot a two-element airfoil so that zero turbulent skin friction is maintained in the 
pressure rise region on the airfoil upper surfaces. Relative movement of the two elements is per- 
mitted m the vertical direction but relative movement in the horizontal direction and deflection are 
not permitted. 

The present paper addresses the problem of determining the position of a leading-edge slat so 
that maximum lift at stall is achieved. It is assumed that the shape of the various elements of the 
lifting system is known. The shapes of these elements are not modified. The theoretical method used 
to calculate the flow field is based on reference 4 and is outlined in Sec. II. Also described in 
Sec. II is the technique used to combine this viscous/potential flow aerodynamic analysis with a 
constrained function minimization code. Section III describes the experimental program which was 
conducted in the Ames 40- by 80-Foot Wind Tunnel to obtain supporting experimental data. The 
comparison of the theory and experiment and the concluding remarks are contained in Secs. IV 
and V, respectively. 


II. THEORETICAL METHOD 


The computerized procedure is a combination of two computer proerams: (1) the aero- 
dynamic analysis program to calculate the flow field and (2) the numerical optimization program to 
determine the slat position. The theoretical basis of these methods is described in the following two 
sections. 


Aerodynamic Analysis 

This analysis provides the aerodynamic information that is required as part of the design pro- 
cedure. In particular, the static pressure distribution on the upper surface of the main airfoil, incipi- 
ent flow separation on the. slat, and incipient flow separation on the main airfoil are used in the 
optimization method and thus must be accurately predicted. The flow-field structure is calculated 
using iteratively coupled solutions of the potential flow and boundary-layer equations. 

The potential flow analysis is an accurate singularity method where the airfoil contours are 
represented by a large number (e.g., 35 to 60 per element) of straight-line seements, on each of 
which is located a linear vorticity distribution. Details of this analysis are presented in reference 4. 
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element. This optimum slat position is determined at .the highest angle of attack at which attached 
How can he maintained on the upper surface of the main element. An alternative view of the analysis 
is that the nearly linear portion of the C\ versus a curve is being extended to as high u lift coeffi- 
cient as possible by proper placement of the leading-edge slat. It is assumed that maximizing the lift 
coefficient at which How separation first occurs will also maximize To . For conficurations that 

stall abruptly, this is a reasonable premise: however, for airfoils that stall gradually, the validity 
o? this assumption is questionable. The validity oi* the method is tested, for a configuration 
with abrupt stalling characteristics, by comparison of theoretical predictions with experimental 
measurements. 

The constrained function minimization method described in detail by Vanderplaats tref. 14) 
and by Vandcrplaats and Moses f-cf. 15) is used in combination with the aerodynamic analysis 
described in the previous section to optimize the slat position. An example of the application of this 
analysis to inviscid single-element airfoil design for minimum wave drag at transonic speeds is pre- 
sented by Hicks et al. (ref. 16). In the present analysis, the function to be minimized is the suction 
peak on the upper surface of the main element. This minimization is performed with the main air- 
foil element at a fixed angle of attack. The design variables are slat horizontal position, slat vertical 
location, and slat deflection. For the three element airfoil considered in this report (see fig. 1 ), the 
slat horizontal and vertical locations are defined by the position of the slat trailing edge. The pivot 
point of the slat is located at the slat trailing edge. Thus, ail three parameters can be varied 
independently. 

Slat translation and deflection are subject to two constraints. The first constraint is that there 
be no tlow separation on the slat upper surface. This prevents the aerodynamic program from 
attempting calculations in the separated flow regime. Mathematically, this constraint is 
Cf s ^ 0.0001, where is the skin friction coefficient at the trailing edge of the upper surface of 
the slat. The second constraint requires that the trailing edge of the slat does not get closer than 
0.02 c to the upper surface of the main element. This prevents the viscous wake of the slat from 
merging with the boundary layer on the wing upper surface. The validity of this constraint is sup- 
ported by the experimental work of Lungstrom (ref. 17) where it is shown that the optimum slat 
gap occurs when the wake of the slat does not interact strongly with the boundary layer on the 
wing. Since the constraint prevents the merging of the slat wake and the boundary' layer on the 
upper surface of the wing, the computationally time-consuming, confluent boundary -layer calcula- 
tion is replaced by a computationally efficient, conventional boundary-layer calculation. This con- 
is expressed as G > 0.02 c, where G is the distance between the trailing edge of the flap and 
the surface of the main airfoil. 

The confluent boundary layer on the upper surface of the flap is analyzed with the slat in its 
initial position. Thereafter, the displacement effect of the confluent boundary layer of the flap is 
held fixed during the numerical search for the slat position that minimizes the suction peak on the 
main airfoil. After the slat position is determined, the confluent boundary' layer is then recomputed 
to ensure that the displacement effect is not substantially different from that with the initial slat 
position. 
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Theoretical Results 


Figure 2 shows a sequence of positions through which the slat moved during the numerical 
design procedure. For this particular example, the angle of attack of the mam airloil is . 
the slat in the initial position (slat position No. 1 ) the slat is lightly loaded, whereas the l«ding-edg 
region of the main airfoil is highly loaded resulting in a large suction peak and a strong adv ersc pres 
sure eradient. After nine iterations, each requiring the calculation of the grad.en of the -uct on 
peak with respect to three slat position variables, and for most iterations the ea culation o e 
Gradient of one or both of the constraints as well, the analysis has converged and the slat posit on 
that minimizes the suction peak on the main airfoil has been determined In this final Position the 
suc'ion peak, and therefore the adverse pressure gradient on the mam airfoil, has been substantially 
reduced in exchange for an increased loading on the slat. The numencal design yields a slat loading 
that is limited only by the constraint that (low separation not occur on the slat upper surface. 
Figure 3 shows the surface pressure distributions for angles ot attack of the main airfoil of 13 , 
] 9 ° an d 24°. At each of these angles of attack, the slat position is such that the suction peak is 
minimized and the boundary layer on the- upper surface of the slat is on the verge of separation. 

Figure 4 shows the variation in slat horizontal position, slat vertical position, and slat deflec- 
tion for angles of attack ranging from 13° to 24°. Also shown is the skin-friction coefficient at the 
trailing edge of the upper surface of the main element. Note that as the angle of attack is increased, 
the boundary layer on the main airfoil approaches separation (Cf m - 0) due to an increase in the 
adverse pressure gradient (Fig. 3) - even though the slat is maintained m a position that minimizes 
the suction peak in the leading^dge region. Extrapolation of these numenca 1 results predicts 
incipient flow separation (Cf m = 0) on the main airfoil at an angle of attack of 25 . The slatposttion 
for this angle of attack is specified as the theoretically predicted optimum location required to give 
maximum lift. This assumption will be verified by comparison with experiment. 


III. WIND-TUNNEL TEST PROGRAM 


As part of this study, wind-tunnel tests were conducted in the Ames 40- by 80-Foot Wind 
Tunnel to verify the reliability of the theoretical predictions by experimentally determining the 
optimum position of the leading-edge slat. 

A rectangular planform wing was used in these tests (fig. 5). The wing was equipped with a 
full-span leading-edge slat and a full-span single-slotted flap. The slat chord was 0.17 c and the flap 
chord was 0.40 c. A detailed description of the slat, main airfoil, and flap shapes is given in refer- 
ence 16. The wing span is 16 m and the extended chord is 2.15 m. The relatively high aspect ratio 
of 7 5 combined with the rectangular planform, resulted in a configuration with neany two- 
dimensional fiow over much of the wing span. Also, the use of a high aspect ratio finite wing elimi- 
nated adverse wind-tunnel wall interference effects associated with two-dimensional, high-lift air- 
foils which span the entire test section. 

The slat brackets permitted continuous adjustment of horizontal location, slat vertical loca- 
tion and slat deflection. The range of movement in the horizontal and vertical directions was 
0.04' c and 0.05 c, respectively. The forces and, moments on the wing were measured along the 
model centerline for a limited number of test conditions. All cata were obtained at a Reynolds 


number of 3.8X I0 fi ami a Mach number of 0.10. The flap location and deflection ancle (5/- 10 ) 
were held fixed throughout the test. The gap between the trailing edge ot the wing and the upp<-r 
surface of the flap was 0.03 c and the Ilap leading edge was located at 0.04 c ahead ot the trailing 
edge of the main airfoil. This flap position is representative of a takeoff flap setting. 

As was discussed in Sec. II. the aerodynamic calculations are based on a method that does not 
model extensive (low-separation effects. Thus, it is expected that this type o! anahsis will be most 
accurate when the configuration being considered has an ,:~rupt. rather than gradual, stall. The mea- 
sured lift curve for this wing is shown in figure 6. It is ne:.:e linear up to an angle ot attack ol -6 . 
At an angle of attack of 26°, the wing abruptly stalls. Thio. flow-separation ellects are expected to 
be negligible at angles of attack below stall. Figure 7 shows a comparison ol the theoretically com- 
puted' and the measured airfoil section pressure distributional a lift coefficient just below stall. The 
experimental data are for a wing angle of attack of 24°. The center-section loading corresponds to a 
section angle of attack of 20° and a section Tg of 3.1. These results verity that the aerodynamic 
analysis method is able to accurately predict the upper surface pressure distributions, even when 
flow separation is incipient. This is essential to the design procedure because at the final design con- 
dition, both the boundary layer on the upper surface of the slat and the boundary layer on the 
upper surface of the wing are on the verge of separating. 


IV. COMPARISON OF THEORY AND EXPERIMENT 


The experimental optimization data are summarized in figures 8-10. The variation of C/, max 
with horizontal j osition of the slat, vertical position of the slat, and deflection of the slat is shown. 
Also shown on the figures are the theoretically optimized trailing-edge position of the slat and the 
slat deflection. The predicted horizontal and vertical locations of the slat are within 0.003 c and 
0.007 c, respectively, of the experimental value. The theoretically determined slat deflection is 
-36°, whereas the experimental value is -40°. The net result is that the theoretically optimized slat 
position is 4° nose up and 0.008 c aft and up relative to the experimentally optimized location. 

Figure 1 1 shows the theoretically and experimentally optimized slat positions relative to the 
leading edge of the main element. With the slat in the theoretically optimized location, the 
measured maximum lift coefficient of the wing is 3.16 and with the slat in the experimentally^ 
optimized location, the measured maximum lift coefficient of the wings is 3.28, a difference ot 
only 4 7c. 

The most significant difference between the predicted and measured optimum positions is the 
4° in deflection. The theoretical results indicate that the slat can sustain a greater loading than was 
found experimentally. Although it is not possible to determine precisely why this occurs, inaccu- 
racy in the prediction of incipient flow separation of the boundary layer on the upper surt ace of 
the slat is probably the cause because this separation has a direct influence on the optimized slat 
deflection. 

The relatively good agreement between theory and experiment shows that the theoretical opti- 
mization technique is useful for estimating the position of leading-edge slats that results in 
maximum lift. Using the theoretical predictions as starting conditions, the wind-tunnel testing 
required to experimentally optimize multielement airfoils can be significantly reduced. 
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V. CONCLUDING RLMARKS 


A design moth- J Tor the optimization of the position of a leading-edge slat so that maximum 
airfoil lift is obtain 'd has been successfully demonstrated. For an airfoil equipped with a leai i ne- 
ed ee slat and a sinele-slotted Hap. the measured maximum lift with the slat in the theorem. alls 
determined optimum position is within 4T of the measured maximum lilt with the slat in the 
experimentallv optimized posi.ion. The analysis is limited to configurations that do not have sig- 
nifieant How ’separation at angles of attack below stall. The details of the How-1, eld structure arc- 
modeled usini! iterativelv coupled potential How and boundary -layer analysis. Accurate prediction 
of incipient flow separation of the upper surface boundary layers on the slat and on the mam wing 
provide the viscous constraints which are essential to the design procedure. It is concluded that, lor 
multielement airfoils that do not have extensive flow-separation ettccts at angles ol attack near stall 
and thus have relatively abrupt stalling characteristics, minimization ol the. suction peak on the 
main component combined with incipient flow separation on both the slat and mam component is a 
valid method for optimizing the position of a leading-edge slat tor maximum lift. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, March 22, 1979 


APPENDIX 


ANALYSIS OF VISCOUS FLOW 


Tllc Levy-Lees form f the incompressible turbulent houndary-layer equations is: 


^ — +l l + /■' = 0 
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The coefficient ? is defined as. 




(A5) 


where 6 is the eddy viscosity. An eddy viscosity 
Cebeci (ref. 1 


model which eombines the recommendation o( 


8) and Cebeci ct’al. (ref. 19) is used to calculate the Reynolds stresses for the trans.- 
•hulent boundary layers. In the inner region of the boundary layer, the eddy viscosity 


tional and turbulent boundary - lay 
is denoted by €('. in the outer regio ^ 

incompressible form and the effects ot mass transfer are neglecte . 
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e outer region it is denoted by e„. When these parameters are reduced to thur 
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where L, the mixing length for the inner region, is given by 
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A parameter y^ r is introduced to control the overall Reynolds stress level during transition. 
At the point of laminar instability, y fr = 0. and for fully turbulent flow, y. r = 1. In the transition 
zone . _ ~ . 

V = 1 ~ CX P (~ (s ~ s ir )G J s (A 1 3) 

where 

(j* Rq 1 - 6 * 

^ = -p- -J- - 5 = 60 + 4.68,V’ (A 1 4) 

s fr and are evaluated at the point of laminar instability. 

The correlations ol Smith (ret. 20) are used to determine the location of laminar instability. 
These equations are 


R fr - -0.4709 + 0. 1 1 060 In Rq - 0.00589 Kin Rq r 

tor Rq < 650. and 


K ir = 0.694 i 2 - 0.23992 In R Q + 0.0205(!n R 0 ) 2 

for Rq > 650 
where 
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indicates laminar instability. 

It laminar or transitional boundary-layer separation occurs, the correlation of Castor (ref. 2 1 ) 
is used tc determine it turbulent reattachment occurs. These equations are 
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sen 


-0.09 


for Rq < 125 , and 


A v<7 , = 0.027 -0.7575X10- 3 R e ~ 1.157X10' 6 Rq 2 


(A 19) 


for Rq7> 125 whore turbulent reattachment is indicated when 

01 ^,-r 

v ds A «v> 


(A20) 


If catastrophic separation is indicated, the boundary -layer calculation is terminated. It laminar sepa- 
ration followed by turbulent reattachment is indicated, then transition to turbulent How is initiated 
at the separation point. Should transitional boundary-layer separation occur, the local adverse pres- 
>ure gradient is reduced until an attached How solution of the boundary -layer equations is obtained. 
In this case, a nearly zero skin-friction solution is obtained for the transitional reattaching boundary 
layer. 
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Figure 1.— Three-element airfoil used in the investigation: 6y= 10°. 







Figure 3.- Airfoil pressure distributions with slat positioned to minimize the suction peak on main 

airfoil: 5,-= 10°. Rc = 3.8XI0 6 . 
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ANGLE OF ATTACK, a, deg 

Figure 4.- Effect of angle of attack on slat position required to minimize suction peak and the 
corresponding skin friction coefficient at the trailing edge of the main airfoil: Sy = 10 , 
Re = 3.8X10 6 . 




Figure 5.— Wind-tunnel model used to experimentally determine optimum slat position required 

for maximum lift. 




LIFT COEFFICIENT, C 



6 _ wing lift coefficient as a function of angle of attack showing nearly linear litt curve 
and abrupt stall: 5 t = -40°, xje = -0.01 ,y s lc = 0.01, 5,= 10°, Re = 3.8X10 6 . 
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Figure 7.- Comparison of theoretical and experimental pressure distributions near stall with a 
section lift coefficient of 3.1: h $ = -42°. x Jc = -0.01, v Jc = 0.02, 5. = 10° 

Re = 3.8X10 6 . s 
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